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ABSTRACT 

We examine the X-ray luminosity of galaxy groups in the CN0C2 survey, at redshifts 0.1 < 
z < 0.6. Previous work examining the gravitational lensing signal of the CN0C2 groups has 
shown that they are likely to be genuine, gravitationally bound objects. Of the 21 groups in the 
field of view of the EPIC-PN camera on XMM-Newton, not one was visible in over 100 ksec 
of observation, even though three of the them have velocity dispersions high enough that they 
would easily be visible if their luminosities scaled with their velocity dispersions in the same way 
as nearby groups' luminosities scale. We consider the possibility that this is due to the reported 
velocity dispersions being erroneously high, and conclude that this is unlikely. We therefore find 
tentative evidence that groups at intermediate redshift are underluminous relative to their local 
cousins. 

Subject headings: cosmology: observations - galaxies: clusters: general - X-rays: galaxies: clusters 



1. Introduction 

According to the paradigm of hierarchical 
structure formation, as matter falls from low- 
density environments to high-density environ- 
ments (i.e., clusters of galaxies) it passes through 
stages of intermediate density, namely, groups. 
Rich clusters of galaxies, containing as many as a 
thousand member-galaxies or more, are visually 
quite prominent and have therefore attracted re- 
search interest for over 70 years. Over the last sev- 
eral decades, however, it has become increasingly 
apparent that small groups of galaxies, contain- 
ing fewer than 50 members and often containing 
as few as 3-10, constitute by far the most com- 
mon environment in which galaxies are found in 
the universe today, and are therefore the domi- 
nant stage of structure-evolution at the present 
epoch. In order to understand the structure and 
evolution of matter in the universe, then, we must 
understand the properties of groups of galaxies. 

In rich clusters, a large fraction of the bary- 
onic mass (> 80 %) exists as diffuse , hot, intra- 
cluster gas - see lEttori et al.l (|2003() and the re- 



view bv lRosati et al.l (j2002l ) and references therein 
- that has been detected with s pace-based X- 
ray telescop es for over 30 years (iGiaccon i et al. 
1974 : Rowan- Robinson fc FabianI 1 1975: Schwartz 
19781 ). If small galaxy groups have a similar hot- 
gas component, then a large fraction of the bary- 
onic ma ss of the universe could be hiding in these 
groups (|Fukugita fc Peebles! I200I and references 
therein). Alternatively, if groups do not contain a 
diffuse, extended intragroup medium (IGM), this 
would be troubling news for the idea that struc- 
ture forms in the hierarchical fashion that is widely 
assumed. 

In the last 15 years, the improved resolution 
and sensitivity of new X-ray telescopes have al- 
lowed us to begin to study the X-ray properties of 
groups of galaxies. ROSAT observations indicate 
that in the local universe, half to three quarters of 



group s h ave a hot. X-ray emittin g IGM (jMulchaev 
|2000( ). IPlionis fc TovmassianI (|2004 hereafter 



PT04) examine the relationship between the X- 
ray luminosity (Lx) and the velocity dispersion 
(a-y) of nearby groups in the lMulchaev et al I (|2003L 
hereafter M03) catalog (although there remain 
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significant uncertainties associated with the anal- 
ysis). At higher redshift, however, a systematic 
study of this kind has not yet been performed, al- 
though some recent evidence suggests that the re- 
lationship at z ^ 0.5 is different from the local one: 
a deep Chandra observation failed to reveal X-ray 
emission from groups discover ed by the DEEP2 
Galaxy Redshift Survey, which Fang et al. ( 2006h 
find strongly suggests that the groups in their sur- 
vey (at z ~ 1) are less X-ray luminous than the 
nearby relationship would predict. 

A recent deep imaging and spectroscopic study 
by the Canadian Network for Observational Cos- 
mology Field Galaxy Redshift Survey {CN0C2) 
made it possible to identify a large number of 
gala xy groups at intermedi ate redshift (0.1 < z < 
0.6) (|Carlberg et al.ll2001allbl ). Gravitational lens- 
ing analysis of the CN0C2 fields yields statisti- 



cal le nsing masses of these groups (jHoekstra et al 
200 if ). This data set, therefore, provides an ideal 



laboratory for studying the X-ray properties of 
groups at these redshifts. 

Fortunately, the data archive for XMM-Newton 
contains approximately 110 ksec worth of observa- 
tions mostly overlapping one of the CN0C2 fields 
(21 of the groups are in the field of view of the 
EPIC~PN camera; see Table[T]). If the relation be- 
tween Lx and ay of groups at these redshifts is the 
same as the relation that holds for nearby groups, 
then we will argue that a few of the most massive 
groups ought to have been visible in the aggregate 
data from these observations. In this analysis, we 
looked for whether the X-ray photons received are 
at all correlated to the groups: (i) Are there more 
(or fewer) photons where there are groups than 
where there are not? (m) Is the spectral energy 
distribution of the photons where there are groups 
different from that of the background? 

In § [51 we describe the data we used and how we 
reduced it. In §[3] and §[11 we describe our analysis 
of the data - spatial analysis and spectral analysis, 
respectively. We found no spatial or spectral evi- 
dence for the groups; and in § [5l we discuss what 
we had expected the data to look like, and how 
surprising it is that we failed to detect the groups. 
In § [SI we attempt to draw conclusions, and we 
speculate as to the reasons why we did not de- 
tect the groups. Finally, in Appendix I, we give a 
brief analysis of several different methods of fitting 
lines to data; and in Appendix II, we describe how 



we estimate the probability that groups at inter- 
mediate redshift share the same X-ray luminosity 
function as those in the local universe. 

2. Observations and Data Reduction 

The CN0C2 survey comprised four fields on 
the sky. We examined the X-ray properties of the 
optic ally identified galaxy gro ups in the 1447+09 
field ( Carlberg et al.l 2001a bl). For our analysis, 
we used the EPIC-PN data from three publicly 
available observations of that field in the XMM- 
Newton data archive, of duration 33 ksec, 33 ksec, 
and 43 ksec. We did not use the corresponding 
EPIC-MOS data in the current analysis because 
of the lower sensitivity of the MOS to low en- 
ergy photons, and the attendant modest increase 
in sensitivity in background-limited images: in- 
cluding MOS would have increased the number of 
counts at the locations of the groups in our anal- 
ysis by a factor of only ~ 50%. The relevant sum- 
mary data on the 21 groups in the field of view 
of EPIC-PN, from lCarlberg et al.l (|2001al ). is pre- 
sented in Table [H 

For basic data reduction, including generat- 
ing good-time-intervals (GTI) files to deal with 
periods of high solar activity, removing known 
hot-pixels, and generating images and exposure 
maps for energy bands of interest, we used XMM- 
Newton Science Analysis Software (SAS) release 
6.1.0. The light-curves used to generate GTI 
files consisted of integrated flux in the range 0.3- 
15.0 keV across the whole detector. Time in- 
tervals of particularly high background (RATE>9) 
were flagged as periods of flaring and were re- 
moved from the data. For the purpose of the anal- 
ysis described below, we considered four energy 
bands: 0.3-0.8 keV, 0.8-1.5 keV, 1.5-4.5 keV, and 
0.2-1.0 keV, which we will hereafter refer to as 
bands A, B, C, and D, respectively. To make im- 
ages in each band, we selected events in the GTI 
with PATTERN<=4 and FLAG==0. Images and cor- 
responding exposure maps were made at 2" per 
pixel, roughly twice the resolution of the intrinsic 
pixel size of EPIC-PN (4.1"). 

Since we were looking for faint emission (from 
extended sources) that was expected to outshine 
the background by only a slim margin, it was 
essential to remove point sources from the im- 
ages. In addition, it was crucial to characterize 



2 



the background carefully before further analysis so 
as to maximize the accuracy of the measurements 
of both the background flux level and the group 
luminosities. Two SAS tasks ~ eboxdetect and 
emldetect - were used to find and remove point 
sources in all observations, in all energy bands; lo- 
cations of point sources were then masked in the 
images and exposure maps. In order to charac- 
terize the background, we used two methods that 
are described in detail in the next section. In 
short, one was to smooth the point-source ex- 
cluded image, and the other was to sample the 
average count-rate over the groups-excluded im- 
age. 

3. Spatial Analysis 

We investigated whether the spatial locations 
of received photons were at all correlated with the 
groups. This inquiry required a metric of surface 
brightness that would allow us to compare quanti- 
tatively different patches of sky. The idea behind 
the desired metric is to find the average number 
of counts per unit sky area per unit time within 
the patch. A reasonable first guess of the metric 
would therefore be the sum of the pixel values in 
the patch, divided by the number of pixels, divided 
by the total exposure time. In an experiment in 
which the pixels have nonuniform effective expo- 
sure time, however, they have unequal sensitiv- 
ity. Maximizing the sensitivity of such an experi- 
ment requires weighting more heavily those pixels 
that are more sensitive in the observation. A rea- 
sonable second guess of the metric of the surface 
brightness of a patch of N pixels, then, is 



11 = kx 



where k is the number of square arcseconds per 
pixel (in our images, fc = 4 arcsec^ pixel"^), pi is 
the number of counts in the ith pixel of the patch, 
Ci is the effective exposure time seen by that pixel, 
and Wi is an arbitrary weighting factor. In order 
to maximize the signal-to-noise ratio of TZ, the 
appropriate weighting factor is Wi = Cj, so in the 
end the metric reduces to 



Tl = kx 



T"" e- 
1^1=1 '^■t 



(1) 



brightness of the groups, we had to know which 
patches of sky to define as the groups. For 
each group, we defined a circular aperture within 
which to search for X-ray emission. The best 
way to assign radii to groups is open to debate, 
so we tried various ways, including two meth- 
ods that depend on the groups' properties - rsoo 
and r2oo, calculated with formulae from the lit- 
erature - and several that do not. The depen- 



dence of rt^nn (lOsmond fc Ponmanll2004l ) and r2oo 
(jMahdavi fc Gellerl 120041 ) on redshift z and line- 
of-sight velocity dispersion are given by 



''soo = 130 kpc 
r2oo = 230 kpc 



100 km s-i 



H{z) 



73 km s-i Mpc" 
Hjz) 

100 km s-iy \73 km s-i Mpc" 



-) (2) 
-) (3) 



where H{z) is the value of the Hubble constant 
at redshift z. All calculations in this paper as- 
sume Ho = 73 km s^^ Mpc~^, ^Im = 0.3, and 
Ha = 0.7. We also tried assigning a fixed phys- 
ical size to each group (125 kpc, 250 kpc, and 
500 kpc). In each case, physical size was converted 
to angular size using a — r/DA{z), where r is the 
group's physical size and Da{z) is its angular di- 
ameter distance. Happily, our main scientific re- 
sults were independent of how we calculated the 
size of the aperture around each group - and this 
is an important point, on which we will elaborate 
below. In this paper, we will present results for 
250 kpc apertures around each group. 

The 7?,~value of a group (according to equa- 
tion ID), within the set of pixels defined by the 
group's angular radius a) cannot by itself be con- 
verted to the group's luminosity, because it in- 
cludes counts from the background that are un- 
related to the groups. The background flux level 
must therefore be calculated and subtracted. The 
true measure of the average surface brightness of 
each group is the difference between the group's 
7?.- value and the 7?.- value of the background: 



SB„ — TZg — TZbg- 



(4) 



In order to determine the average surface 



We calculated the average surface brightness of 
the background in two different ways, and we also 
calculated the local background corresponding to 
each group, which gave us multiple measures of 
the surface brightness SBg of each group. 



3 



Table 1 

The 21 CN0C2 Groups in the Field of View 



Group ID R.A. (J2000) Decl. (J2000) z cr„ (km s^^) Estimated Mass {h~^MQ) 



1 


14 


49 


42 


453 


09 


02 


44 


76 





165 


3 


164 ± 126 


1 


3 


X 


10" 


2 


14 


48 


55 


799 


09 


08 


48 


48 





229 


3 


162 ± 139 


7 


2 


X 


10^2 


3 


14 


48 


57 


854 


08 


57 


05 


58 





262 


4 


229 ± 77 


1 


2 


X 


10" 


4 


14 


49 


28 


406 


08 


51 


54 


07 





270 


3 


104 ± 92 


3 





X 


1012 


5 


14 


49 


13 


637 


08 


50 


10 


34 





271 


4 


112 ± 66 


3 





X 


1012 


6 


14 


49 


44 


220 


08 


57 


37 


84 





273 


3 


165 ± 120 


8 


1 


X 


1012 


7 


14 


48 


49 


049 


08 


57 


50 


15 





306 


3 


93 ± 65 


1 


3 


X 


1012 


8 


14 


50 


20 


387 


09 


06 


00 


62 





306 


4 


199 ± 161 


1 


6 


X 


10" 


9 


14 


50 


13 


703 


08 


57 


37 


28 





325 


4 


175 ± 151 


9 


9 


X 


1012 


10 


14 


49 


03 


618 


09 


07 


02 


86 





359 


4 


82 ± 66 


1 


6 


X 


1012 


11 


14 


49 


40 


860 


09 


02 


15 


68 





372 


3 


126 ± 94 


2 


8 


X 


1012 


12 


14 


50 


22 


850 


09 


01 


14 


74 





373 


4 


44 ± 41 


2 


1 


X 


10" 


13 


14 


49 


49 


201 


08 


51 


23 


08 





374 


4 


291 ± 193 


3 


4 


X 


10" 


14 


14 


50 


00 


825 


08 


49 


06 


65 





394 


4 


308 ± 257 


3 


6 


X 


10" 


15 


14 


49 


14 


155 


09 


11 


15 


38 





394 


3 


394 ± 406 


4 


8 


X 


10" 


16 


14 


49 


55 


433 


08 


56 


11 


39 





394 


3 


507 ± 469 


7 


8 


X 


10" 


17 


14 


49 


32 


746 


09 


03 


41 


28 





468 


4 


488 ± 417 


7 


5 


X 


10" 


18 


14 


49 


29 


974 


09 


09 


08 


15 





469 


4 


217 ± 224 


8 


9 


X 


1012 


19 


14 


49 


31 


428 


09 


05 


05 


22 





472 


3 


123 ± 99 


4 


9 


X 


1012 


20 


14 


49 


23 


516 


08 


58 


47 


83 





511 


3 


565 ± 668 


1 


2 


X 


10" 


21 


14 


49 


23 


763 


08 


55 


17 


84 





543 


3 


151 ± 132 


8 


7 


X 


1012 



''ID Numbers are assigned to the groups that were used for the analysis in this paper. Groups are Usted 
in order of increasing redshift (z). Nz is the number of galaxies with identified redshifts in each group, 
is line-of-sight velocity dispersion of group galaxies . Estim ated mass, and all other numbers in this table 
(except ID Number), are taken from lCarlberg et al.l (2001a). 
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We calculated the average background for the 
whole image in two ways. First, for each energy 
band, we determined the 7?,-value of the image 
for all pixels not contained in any group's aper- 
ture. Second, we created a background image 
for each energy band by smoothing that band's 
image with a Gaussian smoothing radius of 60 
pixels and, as before, determined the 7?.-value 
of the background image for all pixels not con- 
tained in any group's aperture. These two meth- 
ods of determining the average background of the 
whole imag e yielded near ly identical results. Ac- 
cording to IScharj (|2002f ). at redshift of - 0.5, 
low mass objects such as groups should be seen 
most readily in a low-energy band such as D or 
A. Since D-band contains more counts than A- 
band, much of our analysis deals with D-band 
data. The different measurements of the whole- 
detector background were similar; in D-band, for 
example, the background rate was measured to be 
3.26 or 3.08 x 10~^ cts arcsec"^ s~^, depending on 
whether the image or the smoothed (background) 
image was used. 

We also measured the local background cor- 
responding to each group by taking an annu- 
lus of inner radius 1.05a and outer radius 1.20a 
around each group, and finding the 7?.~value of 
the background image in the annulus surrounding 
the group. These local background measurements 
ranged from 2.44 to 4.21 x 10~^ cts arcsec"^ s~^, 



Finally, to understand the spatial variability of 
the background, in each energy band we placed an 
aperture of radius 0.9' (27 pixels) in 20,000 ran- 
dom locations on the point -source-excluded im- 
age, and measured the 7?,-value within each aper- 
ture. This procedure gives an empirical measure 
of the probability distribution function that de- 
scribes the surface-brightness in apertures of ap- 
proximately the sizes of the groups. Results for D- 
band are displayed in the histogram in Figure [T] 
For D-band, the modal 7?.-value of the samples of 
the background is 2.8 x 10~^ cts arcsec"^ s~^; the 
median 7?,-value, 3.0 x 10~^ cts arcsec^^ s^^, is 
similar; and the mean, 3.5xl0~® cts arcsec"^ s~^, 
is somewhat higher owing to a few anomalously 
high measurements that resulted from apertures 
with small numbers of pixels near the edge of the 
detector or near masked regions. 



1200 
1000 



Histogram of /^-Values, 
whole Detector 



10"' 

Value (cts/sq. arcsec/sec) 



Fig. 1. — Histogram of the result of 20,000 random 
samples of the spatial variability of the D-band 
background. The number of samples in each 7?,- 
bin is plotted. The bin-size is 2 x 10~^ counts per 
square arcsecond per second. The mode, median, 
and mean values in this histogram are consistent 
with other measures of the background. 



Surprisingly, the 7?.-values of the groups were 
frequently less than those of the background. 
This means both that 7?,-values of individual 
groups were often less than the correspond- 
ing 7?,-values of the local background, and that 
the stacked 7?,-value for all 21 groups, 2.78 x 
10"^ cts arcsec"^ s""'^, was 10% less than the 
lesser of the two measures of the average back- 
ground. In fact, the average value of the sur- 
face brightness of groups (calculated with equa- 
tion (|4]) and using the local background 7?,- 
value for TZbg] average was weighted by the 
product of the number of pixels and the effec- 
tive exposure time) turned out to be negative: 
—8.6 X 10~® cts arcsec"^ s^^. 

Ultimately, because of the spatial variability of 
the background, which is evident in Figure [l] the 
most useful description of the background is prob- 
ably the local one. The top panel of Figure [2] 
shows the number of counts by which each of the 
21 groups in our sample is brighter than its esti- 
mated local background, plotted against redshift. 
The bottom panel of that figure shows the sur- 
face brightness above the local background SBg, 
again plotted against redshift. The l-cr error bars 
are computed on the ordinates as the summation 
in quadrature of the iV^/^ Poisson noise from the 
group (ag) and from the background (ctbg)- No 
group is seen more than 5ct above the background, 
which is our minimum criterion for detection. The 
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Groups Counts Above 
Background vs. Redshift 



-2c-6 _ 
0.15 



Groups R Above 
Background vs. Redshift 



0.35 
Redshift 



Fig. 2. — Groups intensities above background 
in D-band (0.2-1.0 keV). Top Panel: Counts in 
groups above estimated local background. Bot- 
tom Pannel: Surface brightness of groups above 
background, in units of 7?,- value - i.e., counts per 
square arcsecond per second. None of the groups 
is seen more than 5a above the background, which 
is our minimum criterion for detection. The dot- 
ted line in each panel is the zero line. The dashed 
line in each panel represents the weighted average 
value of the quantity on the ordinate of that panel. 

dashed line in each panel represents the average 
value of the quantity on the ordinate of that panel, 
weighted by the product of the number of pixels 
and the effective exposure time. 

In those cases where both the X-ray and the op- 
tical centers of low luminosity groups are known, 
the correlation betwee n the two shows consider- 
able variance fjMulchaevi |200Q) . Because of the 
possibility that the X-ray centers might not coin- 
cide with the optical centers of the groups in our 
survey, we paid particular attention to the anal- 
ysis with an aperture of radius 500 kpc around 
each group, because an aperture this size should 
be large enough to capture a large fraction of the 
X-rays even if there is an offset of ~ 0.5' between 
the X-ray center and the reported optical center. 
Inasmuch as no group had an 7?.-value Scr or more 
above the local background, the results were sub- 
stantially the same (no plot shown). Moreover, a 
greater fraction of groups had negative S'S-values 
than when the 250 kpc radius apertures were used 
(15/21 vs. 11/21). 

We finally note that the groups might poten- 
tially contaminate the local background estimator. 
Using reasonable models for the intensity and the 



spatial distribution of emission from groups, de- 
scribed in detail in § 15.21 below, we estimate that 
> 30% of flux comes out within our aperture, and 
< 10% of flux comes out in the annulus from 1.05 
to 1.20 times the aperture radius. This portion 
of group flux that falls within the annulus that 
we use to define the local background ought to 
increase the local background measurements, and 
therefore to depress the surface brightness mea- 
surements. The depression, however, is not se- 
vere: since the group surface brightness profile is 
expected to peak toward the center, the surface 
brightness is still expected to be positive. More- 
over, the average 7?.- value of the local background 
samples was 2.9 x 10~^ cts arcsec"^ s~^, which is 
only slightly more than the modal 7?,-value of the 
randomly placed apertures, and slightly less than 
their median value. It therefore does not appear 
that the local background estimator was signifi- 
cantly increased by flux from the groups. 

As a final test of whether the estimate of the 
background was contaminated by group-flux, we 
used an annulus farther from the aperture: the 
annulus from 1.5 to 2.0 times the aperture radius. 
This change had no important effect on our re- 
sults. 

4. Spectral Analysis 

We also explored whether examining the shape 
of the X-ray spectrum would enhance the con- 
trast to the background in the search for groups. 
In the simplest possible test, we investigated 
whether the spectral energy distribution (SED) 
of photons received from the locations of groups 
differed from that of photons from elsewhere. 
The background consists of the particle back- 
ground, and the combined diffuse glows of our 
galaxy, the Galactic halo, the local group, and 
extragalactic point-sources and diffuse sources 
(iMcCammon fc SandersI 1990l : McCammon et al 
20021) . 



Spectral results from our observations are pre- 
sented in Figure [S] In our data, the spectrum of 
the groups is indistinguishable from that of the 
background. This can be seen in several ways. 
The top panel of Figure [3] show, with a bin-size 
of 10 eV, the stacked spectrum of the groups in 
blue, the spectrum of the background in green, 
and the "rest-framed" spectrum of the groups in 
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Groups 

Background 

Groups. "Rest-Frained" 




Fig. 3. — Spectra showing 7?.-values per 10 eV 
bin. . The rest -framed groups show no signs of any 
spectral features Top Panel: Groups, background, 
rest-framed groups spectra, 0.2-2 keV. Spectrum 
for locations of groups is in blue, for background 
is in green, and for rest-framed groups is in red. 
Bottom Panel: Difference spectrum, 0.2-2 keV. 1- 
a error bars are in green. 

red (that is, the energy of each photon from the 
groups is multiplied by (1 + Zgioup) prior to bin- 
ning). The groups and the background both show 
broad emission features centered at ~ 0.55 keV 
and ^ 1.5 keV (top panel of Figure [3]), that are 
most likely due to O Vll at redshift and neutral 
Si on the detector, respectively. The difference 
spectrum, however, is to within the l-cr error 
bars (bottom panel) , indicating that the SED from 
regions of groups is identical to that of other re- 
gions. Furthermore, the rest-framed groups spec- 
trum shows no signs of any spectral features. 

5. Discussion 

In recent years, research efforts have increas- 
ingly been devoted to characterizing the extended 
X-ray emission from collapsed halos of a vari- 
ety of masses, from clusters to groups to galaxy- 
mass halos. For two examples (among many) of 
studies of gas in h igher- mass systems - clusters, 
Akahori fc Masai ( 2005f) investigate the Lx~-Tx 
scaling relation describing intra cl uster gas, and 
Arnaud. Pointecouteau. &: PrattI (j2005l ) investi- 
gate the M-T scaling relation of intracluster gas. 

Our study is concerned wit h the prope r ties o f 
gas in lower-mass systems. Xue fc Wu ( 2000[ ) 



attempt to bridge the gap between clusters and 
groups in the local universe by performing an anal- 
ysis of the Lx-T, Lx-Ov, and a^-T relations of 66 
groups and 274 clusters taken from the literature. 
The XMM-Large Scale Structure survey identified 
gro ups and clusters out to redshift 0.6 and beyond, 
and Willis et al.l (2005) examine the Lx-Tx rela- 
tion of the objects in their survey. Interestingly, 
in contrast to the analysis in this paper, they de- 
tect X-ray emission from low-mass objects out to 
z = 0.558. The groups that they found were X-ray 
selected, so the selection effect surely contributed 
to their finding groups in their survey while we did 
not in ours (although we note that this does not 
explain the physical difference that is responsible 
for the discrepancy in X-ray luminosities). 

5.1. Fits to Lx—CTy Data 

In order to estimate how our data ought to have 
appeared, we need to know what X-ray luminos- 
ity should be expected from a group with a given 
mass or velocity dispersion. The groups of the 
MOB catalog are a convenient sample of nearby 
groups with known velocity dispersion and X-ray 
luminosity. Since any fit to the Lx-(Jv relation in 
the M03 groups will show scatter, a fit is not a line 
or curve but a 2-dimensional region in Lx versus 
ay space. 

In fact, properly, one would want to construct 
a 2D probability distribution function on Lx-cTy 
space based on the distribution of the M03 sample. 
This function would give the probability density 
for a group drawn from the same population as 
the reference sample to have a particular velocity 
dispersion and X-ray luminosity. One could then 
properly examine the likelihood that the CN0C2 
groups are drawn from the same population as the 
MOB groups. 

In practice, it would be difficult to construct the 
2-D distribution function we require, because the 
M03 sample is not dense enough in Lx-cTy space. 
We therefore reverted to a conventional paramet- 
ric analysis, based on linear fits to (logcTt,, \ogLx) 
data. We will present several linear fits to the 
M03 data, and we describe in detail our choice for 
which linear fit is most appropriate. 

PT04 examine the relationship between the X- 
ray luminosity (Lx) and the velocity dispersion 
(cTt,) of galaxy groups in the M03 catalog. They 
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perform regression analysis on logLx and logCTi, 
and find that, as with any regression in which 
there is significant scatter, which hne is considered 
the best-fit depends on which variable is consid- 
ered to be independent and which dependent. We 
repeated their analysis of the M03 groups and ob- 
tained best-fit lines that are very similar (though 
not identical) to those that they obtained. In this 
paper, we present our own fits to the data. If ve- 
locity dispersion is taken to be the independent 
variable, the best-fit regression line is 

log(Ljf/l erg s-i) = 36.8+2.08 log (ct„/1 km s'^) 

(5) 

If, on the other hand, luminosity is taken to be 
the independent variable, the best-fit regression 
line (what PT04 refer to as the "inverse regression 
line") is 

log(Lx/l erg s"^) ^ 30.9-^4.50 log (ct^/1 km s^^) 

(6) 

Data from M03 and best -fit lines are displayed in 
Figure U) 

Which regression represents the "true" relation 
between the X-ray luminosity and the velocity dis- 
persion of groups of galaxies? In Appendix I, we 
present a more detailed analysis of the relative 
merits of various procedures for fitting lines to 
data. The result of this analysis is that neither di- 
rect nor inverse regression is ideal, because surely 
the measurements of both and Lx were subject 
to errors, and furthermore there is undoubtedly in- 
trinsic scatter in the Lx-Uy relation. The better 
relation out of these two is the one in which the 
variable taken to be independent has lower rela- 
tive uncertainty. In our case, this is probably the 
inverse regression, because a velocity dispersion 
determined from only a few galaxies is likely to 
be quite uncertain. Furthermore, the steeper (in- 
verse) relation comes closer to fitting the high-cr^,, 
high-Ljf end of the distribution of groups, which 
is the end that we are most interested in when 
looking for groups at z ~ 0.5. Still, although the 
inverse regression is probably preferable to the di- 
rect regression, since there is scatter in both vari- 
ables, this situation calls for least squares orthog- 
onal distance fitting (described in Appendix I). 
If we assume that the typical relative errors on 
log((T^) and log{Lx) are the same (i.e., the error 
is shared equally between the two variables), then 



the distance-fit line, 

log(Ljf/l erg s-^) = 31.5+4.25 log (ct„/1 km s~^) 

(7) 

is nearly as steep as the inverse regression line, as 
shown in Figured 

5.2. Limits on Group Luminosities 

Since none of the groups was visible at 5a above 
the background, we may place limits on their 
fluxes and therefore, since we know their redshifts, 
on their luminosities. The total number of back- 
ground counts TCbg the location of a group is 
estimated from the smoothed background image. 
Assuming Poisson statistics, the standard devia- 
tion on this number is otCbg ~ {TCbg)^^^ ■ The 
maximum number of group counts TCq, then, is 
limited by 



TCg < ^gtCbg = ^vTCbg^ 



(8) 



We estimate each group's tem perature with th e 
standard relation given in, e.g., iMulchaevI (|2000f l: 



(9) 



where nip is the proton mass, /it is the mean molec- 
ular weight, k is Boltzmann's constant, and /?, 
whose empirical value tends to be around 2/3, 
is the eponymous parameter of the isothermal-/? 
model for density and may be considered to be de- 
fined by this equation. (For low-density objects, 
/3 has been measured to be closer to 1, but for 
analytical simplicity we used (3 = 2/3 for most 
calculations.) The /3 density model for a spherical 
isothermal plasma is 



p(r) 




-3/3/2 



(10) 



where po is the central density, r is the radial 
distance from the sphere's center, and r,. is the 
"core radius," a distance within which the den- 
sity is nearly constant (approximately po). Using 
the Astrophysical Plasma Emission Code (APEC) 
model in XSPEC, we solve for what X-ray lumi- 
nosity, at a given temperature, would lead to the 
number of group counts given by the maximum 
possible value of TCg in © ■ 

In Figure [31 the X's show the maximum lu- 
minosity that each group could have had during 
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our observations (if the group were more lumi- 
nous, it would have shown up at > 5a above the 
background). If a group is predicted to be less 
luminous than the luminosity represented by its 
X, then it is no surprise that we failed to see the 
group; conversely, if a group is predicted to be 
more luminous than the luminosity represented by 
its X, then the group is less luminous than pre- 
dicted. The dashed line shows the Lx-dy direct 
regression line, the dashed-dotted line shows the 
inverse regression line, and the solid line shows 
the distance-fit line. All X's (all groups) lie above 
the direct regression line, so if that line represents 
the true relationship between the variables, there 
is no surprise in seeing none of the groups. Three 
groups, however, lie below both the inverse regres- 
sion line and the distance-fit line; so if those lines 
represent the true relationship between the vari- 
ables, it might be surprising that we failed to see 
three of the groups. In the remainder of this sec- 
tion, we will try to quantify how surprising it is. 

5.3. Quantifying Surprise 



The CN0C2 catalog of small groups (jCarlberg et al 



2001bl ) contains groups with as few as three galax- 
ies in them. Any measurement of the velocity 
dispersion of such a group will necessarily have 
large uncertainties. In fact, for some groups, the 
estimated uncertainty in velocity dispersion is ac- 
tually greater than the velocity dispersion itself. 
Because the symbol a is frequently used to denote 
both uncertainty and velocity dispersion, it is in- 
elegant to represent the uncertainty in velocity 
dispersion, but, eschewing elegance, we shall use 
(Tcr„ to denote uncertainty on ay. 

It is possible that our null result - our fail- 
ure to detect any of the groups in our survey 
- occurred because the velocity dispersions of 
the three groups predicted to be observed were 
overestimated. If so, then perhaps their pre- 
dicted luminosities based on, e.g., ([7]), should have 
been low enough that the groups should not have 
been seen after all. Of course, if there were no 
bias to the measurements of ay, it is just as 
likely that the velocity dispersions of these three 
groups were underestimated, which would make 
their non-detection even more surprising. Fur- 
thermore, in order for this sort of statistical error 
to explain our null result, not only would the ve- 
locity dispersions of the three groups that we think 
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Linear Fits to M03 Group-Data 
And Empirical Upper Limits On 
CNOC2 Group-Luminosities ° 



Direct 

- - Inverse 

L.S.O.F 

X E.U.L. 
o M03 



CT^ (kni/s) 

Fig. 4.— Linear Fits to M03 Group-Data 
And Empirical Upper Limits On CN0C2 Group- 
Luminosities: bolometric X-ray luminosity vs. 
line-of-sight velocity dispersion. The circles show 
the M03 groups. The dashed line is the direct re- 
gression line, the dashed-dotted line is the inverse 
regression line, and the solid line (L.S.O.F.) is the 
least squares orthogonal fit line. These three lines 
are fits to the M03 data. The X's (E.U.L.) dis- 
play the "empirical upper limits" to groups' lumi- 
nosities, i.e., the maximum luminosities that the 
CN0C2 groups could have had during our obser- 
vations without showing up at more than 5a above 
the background. The uncertainty in velocity dis- 
persion associated with the X's is typically nearly 
~ 100%, as shown in Table [TJ Three of the X's 
lie beneath both the dashed-dotted and the solid 
lines, indicating that these three groups were pre- 
dicted to be luminous enough to be detected, and 
yet were not. 

we ought to have detected need to have been suf- 
ficiently overestimated, but no other groups could 
have have had their velocity dispersions underesti- 
mated by too much. In short, the statistical errors 
would need to be in particular directions on par- 
ticular groups. If there is no non-statistical reason 
for measurements to be in error, then what is the 
likelihood that the types of error needed to explain 
our null result actually occurred? 

To answer this question, we simulated 10,000 
mock-CiVOC2 catalogs, each with its own set 
of mock-velocity dispersions {{a^{i)}) associated 
with the mock-groups. If the set of measured ve- 
locity dispersions from the actual CN0C2 catalog 
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is {ay{i)}, then, for each mock-catalog, the sim- 
ulated velocity dispersion of the ith group 0^(1) 
was drawn from from a gamma distribution with 
mean av{i) and standard deviation a^^ii)- The 
gamma distribution has the attractive property 
that the mean and the standard deviation may 
be independently specified. Because the tail of 
the gamma distribution extends to infinity, how- 
ever, a small fraction of mock-groups were initially 
given unrealistically high velocity dispersions. In- 
cidentally, this is preferable to the normal distri- 
bution, which has another tail that extends below 
zero to negative infinity - unphysical values for 
velocity dispersion. To deal with the gamma dis- 
tribution's problem of predicting a few unrealisti- 
cally high values, we, for each group i, eliminated 
from the set of mock-catalogs all mock-groups 
with velocity dispersions more than 4crcr^ (i) above 
the mean (ct„(i)), and replaced these values with 
new gamma-distributed variables. We repeated 
this procedure iteratively until there were no more 
mock-groups with velocity dispersion more than 
4crcr„ above the mean. Because this procedure in- 
volves replacing the high-end outliers with values 
closer to the mean, it had the effect of slightly 
depressing the mean. All values were therefore in- 
creased by the difference between the new mean 
and ay{i), thereby raising the mean back to cru(?), 
as it should be. 

We also tried several other random distribu- 
tions for assigning mock-velocity dispersions to 
our mock groups. We tried a gamma distribu- 
tion with no high-end cut-off; and we tried sev- 
eral modifications of the normal distribution, all of 
which included a low-end cut-off at 0. We found 
that those results were not importantly different 
from results from the gamma distribution with the 
high-end cut-off, which is the analytically sim- 
plest distribution that produces realistic results. 
Results are presented only for the mock-catalogs 
generated with the high-end cut-off gamma dis- 
tribution. 

Each mock-group was assigned the same aper- 
ture size as the corresponding real group. Its 
mock-intensity was spatially distributed accord- 
ing to a spherical isothermal-/? model. Since emis- 
sivity is proportional to the square of density p, the 
surface brightness or intensity / is proportional to 
the integral along the line-of-sight of p^. At a 
projected viewing "impact parameter" b from the 



sphere's center. 



I{b) cx 



pirfds, 



where r^, the "halo radius," is the outer bound- 
ary of the sphere of plasma, and s is the variable 
along the line-of-sight. Substituting for p with 
equation (fTU)l and performing the integration re- 
sults in the following expression for the surface 
brightness: 

(11) 

where x = b/rc and A = rh/vc are substitutions 
to simplify the integration, and 2-^1 is the Gauss 
hypergeometric function. Three quantities go into 
the normalization of B for each mock-group: its 
mock-luminosity from e.g. ([7]); its luminosity dis- 
tance Dl{z); and its core radius. Finally, we 
assumed a foreground colum n density of neutral 
hydro gen of 2 x 10^° cm^^ (jPickev fc LockmanI 

It was unclear what core radii our mock-groups 
should be assigned, because the /3-model core 
radii of groups and clusters vary a great deal. 
Ota fc Mitsudal (j2004l ) report that in a sample of 
79 clusters, including 45 "regular" clusters and 
34 "irregular" clusters, the mean of all 79 val- 
ues of Tc is (163 ± 202) /i7o~^ kpc, midway be- 
tween the mean Vc of the regular clusters ((76 ± 
60) /i7o~^ kpc) and the mean Vc of the irregular 
clusters ((273 ± 259) hro~^ kpc). Because there 
is no strong correlation between Vc and group- 
mass, we think it is reasonable to adopt a uni- 
form assumed value of Vc for the mock-groups in 
our mock-catalogs. We ultimately decided to set 
each mock-group's core radius to 125 kpc, because 
this value lies i n the middle of t he ran ge of val- 
ues reported bv lOta fc Mitsudal (|2004l ). So long 
as the value of fc that we adopted is physically 
plausible, the precise value does not matter much, 
because the predicted number of counts from a 
group, within an aperture of fixed size, is not very 
sensitive to the assumed value of Tc, as long as the 
angular-size of the aperture is larger than that of 
the core radius. 

The SED of each mock-group (and hence the 
flux in a given energy band) is calculated with 
XSPEC, using APEC, with its mock-temperature 
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Fig. 5. — Histogram of the results of the mock- 
catalog analysis. Of the 10,000 mock-catalogs, 
96.5% contained at least 1 mock-group visible at 
the 5a level. In the actual CN0C2 survey no 
groups were visible at this level. 

calculated from equation For A, we tried val- 
ues between 2 and oo, and for /3, we tried values 
between 0.6 and 1.0. The results presented be- 
low are for A = oo and (3 = 2/3, and for D-band 
(0.2-1.0 keV). 

In brief, the results of this Monte Carlo indi- 
cate that the probability appears to be small, but 
not vanishingly so, that statistical fluctuations in 
the measured values of ay explain our null result. 
Figure [5] shows a summary of the mock-catalog re- 
sults. The median and the mode number of groups 
seen per catalog are both 3, and the mean number 
seen is 3.1, which is in good agreement with the 
prediction of section 15.21 The finding of our ob- 
servation - that no groups were seen - occurred in 
only 3.5% of the mock-catalogs. If there were no 
systematic errors in the measurements of velocity 
dispersions of the CN0C2 groups, then it is un- 
likely but not impossible that the only reason we 
had expected some of the groups to be visible in 
the first place was erroneously high measurements 
of ay. 



5.4. Redshifted Mulchaey Groups 

Perhaps wc have placed too much emphasis, in 
sections 15.11 and 15. 2i on which best-fit line actu- 
ally fits the M03 data the best - when the truth 
is that there will be significant scatter to the data 



around any line drawn through the cloud of points. 
An alternative way to determine whether it is sur- 
prising that no groups in our survey showed up 
at 5(7 above the background is to ask whether the 
CN0C2 groups and the M03 groups were drawn 
from the same population. We addressed this 
question by investigating how many of the groups 
in the M03 atlas would have been visible at the 5a 
level if they had been in our observation. 

The M03 groups are all in the local universe, 
approximately at redshift 0. If they had been in 
our survey, i.e., between rcdshifts 0.1 and 0.6, they 
would have both been dimmer and subtended a 
smaller solid angle on the sky. Furthermore, they 
would have been observed with XMM~Newton, 
which, owing to its greater focal ratio, has a higher 
particle background rate than ROSAT (the instru- 
ment with which they were originally observed). 

There were 109 groups in the M03 catalog, of 
which 61 were seen in X-rays. We simulated ob- 
servations of each of these 61 M03 groups at each 
of the 21 redshifts of the CN0C2 groups, rang- 
ing from zi = 0.165 to Z21 = 0.543 (ID numbers 
taken from Table [T]). As before, our criterion was 
that a group must outshine the background by at 
least 5a to qualify as being detected in our sim- 
ulated observations. Each group was assumed to 
be at a location with detector-averaged charac- 
teristics: the mean exposure map value, the me- 
dian background value (from Figure [5]) , the mean 
amount of smearing due to the point-spread func- 
tion, and the mean amount of lost usable detector- 
area because of bright point sources. Much like in 
the mock-catalog calculations described in section 
5.31 we calculated the count-rate of each group 
with XSPEC, using APEC, with the group's tem- 
perature calculated from equation ([9]). We dis- 
tributed emission according to the isothermal-/? 
profile using equation (fTT|) . with f3 = 2/3, A = 00, 
and Vc — 125 kpc. For our simulated observa- 
tions, we used an aperture of 250 kpc, and looked 
in D-band. For background, we used TZbg = 
3.0 X 10"^ cts arcsec"^ s~^. Finally, as before, we 
assumed a foreground column density of neutral 
hydrogen of 2 x 10^° cm~^. Within our simula- 
tion, the number of M03 groups that exceeded the 
5a detection criterion varied from 22 of the 61 to- 
tal groups (or pi = 36%) at the lowest redshift 
(zi), to 2 of 61 {p2i = 3%) at the highest redshift 
(^21). 
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We identify two populations of the Mulchaey 
groups - all 109, and the 61 that were seen in X- 
rays. In Appendix II, we present in detail our esti- 
mate of the likelihood that the CN0C2 groups and 
the M03 groups were drawn from the same pop- 
ulation. The result is that the probability of ob- 
taining our null result on the assumption that the 
CN0C2 groups were drawn from the same popu- 
lation as the whole M03 catalog is 14.7%; and the 
analogous probability on the assumption that the 
CN0C2 groups were drawn from the same popu- 
lation as the 61 X-ray-bright M03 groups is 2.8%. 
These probabilities (especially the latter) are low 
enough to be interesting, but are certainly not dis- 
positive. 

6. Conclusion 

The analysis in section [53] is mildly suggestive 
that the groups of the CN0C2 survey constitute 
a sample that is on average less luminous than 
the sample of groups observed by Mulchaey and 
collaborators in the 2003 atlas. The interesting 
comparison, however, is not between two particu- 
lar samples, but rather between the population of 
groups at intermediate redshift and the population 
in the local universe. 

A first caveat is that the two samples were se- 
lected in very different ways, and there is no a pri- 
ori reason to suspect that their underlying popula- 
tions have similar luminosity functions. Further- 
more, it is conceivable that for some reason the 
M03 groups should be expected to be a more lumi- 
nous sample of their population than the CN0C2 
groups are of theirs. Perhaps the M03 groups 
are on average more massive (although the ranges 
and distributions of velocity dispersions for the 
two samples are nearly identical). Do the M03 
groups contain more galaxies than the CN0C2 
groups? Certainly, Mulchaey reports more galax- 
ies per group on average - including as many as 
63 galaxies in one group. Of the CN0C2 groups 
in our field of view, none contained more than 4 
galaxies with identified redshifts. But most or all 
of the difference here is almost certainly due to 
the difficulty of measuring spectroscopic redshifts 
of faint galaxies at large distances; we therefore 
see no reason to think that the M03 groups are on 
average more massive than the CN0C2 groups. 

An obvious criticism of the present work is that 



the CN0C2 groups might not be groups at all, but 
just chance superpositions of galaxies. While this 
would clearly explain the lack of X-ray emission, 
we will argue that we find it unlikely to be the 
reason why we detected no hot gas. 

ICarlberg et aD ()2001bl) identified as a group 
any object in the CN0C2 Survey that contains 
at least 3 galaxies in sufficiently close redshift- 
space proximity. This identification procedure 
does leave open the possibility of confusion: if 
two galaxies are separated by a large distance, but 
have appropriately large peculiar velocities, they 
could have similar recessional velocities (redshifts) 
despite not being part of a gravitationally bound 
structure. Or, even if a few galaxies are close to 
one another in physical space, they will not be part 
of a bound system without a dark matter halo that 
is massive enough to bind them. 

In searching for groups and clusters, it is of- 
ten assumed that where there are associations of 
galaxies, there is dark matter to bind them, but 
this assumption could be wrong. With rich clus- 
ters, there are multiple independent checks of this 
assumption: the temperature of intergalactic gas 
(which can be detected through its X-ray emission 
or through inverse Compton distortions to the cos- 
mic microwave background) is one indicator of the 
mass; the velocity-distribution of galaxies is an- 
other indicator of the mass, and in particular the 
distribution should be Maxwellian (or Gaussian) 
for a virialized structure; and the degree of grav- 
itational lensing of background galaxies is a third 
indicator of the mass. 

The first check has obviously failed in this case 
- for these 21 groups, there is no apparent hot 
gas component. The second check is useless for 
an association of very few (3-6) galaxies - nearly 
any velocity distribution is consistent with an un- 
derlying Gaussian distribution. How about lens- 
ing? For structures of ^ IO^^'^^Mq, noise prevents 
detection of a lensing signal. As a result, in the 
CN0C2 fields, it is only possible to determine a 
statistical lensing mass for the groups by stack- 
ing them and determining their combined lens- 
ing signal. Th i s analysis h as been performed by 



Hoekstra et al.l (|200lL 120031 ): their calculated esti- 



mate of the average velocity dispersion per group 
based on the mass inferred from the average shear- 
ing of background galaxies (258j^5g km s~^ for 
= 0.2 and Ha = 0.8) is in good agreement with 
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the average velocity dispersion of 230 km s~^ from 
the kinematics of the groups' member-galaxies. 
The observed lensing signal was sufhciently strong 
that chance projections of unbound galaxies alone 
would have difficulty mimicking it. Some of the 
detected lensing mass is therefore almost surely 
in the form of galaxies orbiting within a parent 
dark matter halo - groups. It is entirely reason- 
able to expect part of this group-component of the 
lensing mass to consist of hot gas trapped within 
the group potential. On the other hand, it is not 
unlikely that a portion of the lensing mass is in 
the form of projected, yet physically unassociated, 
galaxies. Although there may be limited amounts 
of gas associated with the individual galaxies, such 
projections will produce no extended group emis- 
sion. Therefore, the presence of a lensing signal 
does not automatically lead to the expectation 
of group-like X-ray emission that is at sufficient 
surface-brightness to be detectable in our survey. 

It is conceivable that there were several rel- 
atively massive groups, or dark clusters (i.e., 
cluster-mass objects with very low mass-to-light 
ratios), that dominated the mean mass measure- 
ment. Although we cannot rule out the possibility 
that a large fraction of the mass is in only a few 
objects that happened to fall outside the field of 
view of our observations, and that the remaining 
objects are not gravitationally bound structures 
but rather chance associations of galaxies that 
just happen to be near one another in redshift- 
space, this seems somewhat unlikely: Combinato- 
rial statistics dictate that if there were a few mas- 
sive halos scattered among the 40 alleged groups 
and none of the rest are true groups, then the 
probability is low that all the massive ones ended 
up among the 19 groups that did not fit into our 
field of view. Consider, for instance, if there were 
only four massive halos: since 19 is roughly half 
of 40, if the four massive halos were randomly 
distributed the probability that all four would 
be out of our field of vi ew is ^ (1/2) '^ Fur- 
thermore, a calculation in lCarlberg et al. ( 2001bl ) 
shows that the number of alleged groups in the 
field in the redshift range 0.1 < z < 0.6 (40 in the 
CN0C2 field, 21 in the EPIC-PN field of view) 
is roughly consistent with predictions. The con- 
cordances both of the observed number of groups 
with the predicted number, and of the observed 
velocity dispersions with the mass inferred from 



weak lensing, together constitute a strong indi- 
cation that the CN0C2 group sample does not 
suffer significant contamination from chance su- 
perpositions of galaxies. An important caveat is 
that the number of groups in the CN0C2 fields 
in the high range is in excess of the aforemen- 
tioned predictions derived from t he basic formal- 
ism of IPress fc S chechtei l (Il974h W a factor of 
~ 4. Based on the Carlb erg et al.l (j2001al ) results, 
~ 80% of such high (t„ groups could be false, and 
our own null results are consistent with such this 
hypothesis. 

If we take seriously the conclusion, for which 
our analysis has generated some modest evidence, 
that galaxy groups at intermediate redshift are less 
X-ray luminous than those at redshift zero, we can 
think of several potential explanations. Groups at 
z ~ 0.5 are younger than their relatives in the local 
universe, which leads to the following three pos- 
sibilities: (i) Since simulations and observations 
indicate that groups and clusters at the present 
are accreting gas from cosmic filaments, groups at 
a younger stage of the universe would presumably 
have been accreting gas for less time than the more 
evolved groups of the present-day universe and 
might therefore contain less gas. (ii) Even if the 
intragroup gas in young groups is not lower density 
than the gas in groups today, if the young groups 
are still in the process of heating up - if they have 
not yet reached their virial temperatures - their lu- 
minosities would be reduced. (We find this possi- 
bility to be unlikely, because the heating timescale 
is short relative to a Hubble time, even at z ~ 0.5.) 
(iii) Even if the gas in young groups is on aver- 
age the same density and the same temperature as 
the gas in nearby groups, if it is of lower metallic- 
ity its X-ray emissivity would be reduced. There 
are undoubtedly feedback processes that enrich 
the intragroup medium, and pristine gas of pri- 
mordial abundance is also surely constantly falling 
onto groups, so it is not clear whether to expect 
higher metallicity in the IGM of groups that are 
at z ^ 0.5 or at z ~ 0. For warm gas, observed 
in a soft X-ray band such as our D-band, line- 
emission can be important, depending on abun- 
dances. For example, if a T = 1 keV group at 
redshift 0.5 contained primordial gas, its D-band 
X-ray flux would be only two thirds as great as if 
its gas were at 10% solar abundance and only half 
great as if its gas were at 20% solar abundance. 
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It is possible, therefore, that the relation between 
velocity dispersion and X-ray luminosity ought to 
be shifted down by a factor of ~ 2 in Figure [H 
Such a shift would make it less surprising that we 
detected no X-ray emission from the groups. 

Future data on high-redshift clusters that make 
it possible to examine the time-evolution of the 
density and metallicity of the intracluster medium 
in low-mass clusters, and deeper searches for X- 
ray gas in low-mass objects at intermediate red- 
shift, will help to answer why our survey failed to 
detect such gas. 
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Appendix I 

In this Appendix, we address the question of 
which regression represents the "true" relation be- 
tween the velocity dispersion and the X-ray lumi- 
nosity of groups of galaxies. It is instructive first 
to consider a statistical effect that is common to 
all regression analysis. It has been known at least 
since Karl Pearson's seminal paper in 1901 that 
if there is scatter (say, due to measurement error) 
in the variable considered to be the independent 
one, this will in general cause the recovered slope 
from regression analysis to be shallower than the 
slope of the true underlying relation between the 
variables. (PT04 demonstrate a nice example of 
this effect.) An important consequence of this ef- 
fect is that when there is scatter in both variables, 
the slopes from both the direct and the inverse re- 
gression will be shallower than the respective true 
slopes. 

For example, suppose two variables x and y are 
related hy y = ax (or x = /3y, where /3 = 1/a). 



In an experiment, a large number of (x, y) pairs 
are measured, and there are statistical uncertain- 
ties in the measurements of both variables. If x is 
considered to be the independent variable, regres- 
sion analysis will lead to y = ax, where in general 
a < a because of the scatter in x. If y is con- 
sidered to be the independent variable, regression 
("inverse regression") analysis will lead to x = f3y, 
where in general (3 < [3 because of the scatter in y. 
Of course, if we wish to represent y as a function of 
X even in the case of the inverse regression, we will 
write y = ax, where a = 1/ (3 > a. So, although 
the recovered slopes are in general less than the 
respective true slopes for both regressions {a < a 
and (3 < (3), when we fix the abscissa and ordi- 
nate variables we get a shallower slope than the 
true one in the case of the direct regression and 
and a steeper slope in the case of the inverse re- 
gression. Note that this effect holds whether the 
scatter comes from measurement error or is intrin- 
sic to the variables. 

When there is scatter in both variables, so that 
neither one is a proper independent variable, it 
is sometimes appropriate to use a type of best-fit 
line that depends on the geometric orientation of 
the points. Instead of finding the line that mini- 
mizes the summed, squared vertical deviations of 
the data from the line, as in the case of regres- 
sion, some situations are better analyzed by find- 
ing the line that minimizes the summed, squared 
perpendicular distances of the data from the line. 
This is sometimes called "least squares or thogonal 
distan ce fitting," and was the subject of !Pearso3 
(jlQOll ) . Because this type of fitting, in contrast to 
linear regression, is unfortunately sensitive to the 
units of measurement, it is advisable to normal- 
ize the data prior to fitting. The most common 
normalization is by the estimated erro rs or uncer- 
taint i es (cr.r and (7„) on the varia bles (jPress et al 
19921 : lAkritas fc Bershadvl [l996h . The direct re- 



gression line will have the minimum slope, the in- 
verse regression line will have the greatest slope, 
and the "distance-fit" line will have an interme- 
diate slope (closer to the direct regression slope if 
the scatter in x is smaller; closer to the inverse re- 
gression slope otherwise). All three lines intersect 
at the centroid of the data. 

For completeness, we present the equation of 
the least squares orthogonal distance fit line. Let 
us express the measured data as a set of N points 
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{{xi,yi)}^i. To simplify the form of the hne's 
equation, we define the following sums of data: 



N 



N 



Sx — Xi Sy — ; 



N 



N 



N 



^ ^ Xi Sxy — ^ ^ (•^i Hi) ^y^ — ^ ^ 



1=1 



i=l 



(12) 

and we define the following combinations of these 
sums: 

A = SxSy NSxy 

B = [Sxf-{Syf-N{Sx2^Sy2) 

C = -SxSy+NSxy. (13) 

If the equation of the least squares orthogonal dis- 
tance fit line is 

y = mx + 6, (14) 
then the slope of the line is given by 
-B ± - 4AC 



2A 

and the y-intercept is given by 

Sy VflSx 

6 = 



N 



(15) 



(16) 



The expression for m in equation (|15p contains a 
±-sign from the solution to a quadratic equation. 
One solution for m gives the line that minimizes 
the summed squared orthogonal distances to the 
points, and the other solution gives the line pass- 
ing through the centroid of the data that maxi- 
mizes the summed squared orthogonal distances 
to the points. 

Appendix II 

In this Appendix, we address how to compute 
the probability that the CN0C2 groups and the 
Mulchaey (M03) groups derive from the same pop- 
ulation. In order to do this, we find the proper- 
ties of the hypothetical parent group population V 
that maximize the joint probability of our results 
and Mulchaey's results. 

Let the redshifts of the CN0C2 groups be la- 
beled {zk}l^i- At each redshift Zk, let Nk be the 



number of M03 groups that our simulation pre- 
dicts would have been visible in our 110 ksec obser- 
vation, and let pk be the probability that a group 
drawn from V will be visible in our survey if it is 
in our field of view and at redshift z^. 

At redshift Zk, the probability of obtaining our 
result - of not seeing any groups - is 



qk = (1 -Pk); 



(17) 



and the probability of obtaining the M03 result - 
of seeing Nk groups, as predicted by our simula- 
tion - is 



B{M,N,,p,)^f^^)p,^^qk''-''\ 



(18) 



where M is the total number of groups in the M03 
sample, either 109 or 61, as described above. 

The joint probability, then, of obtaining our re- 
sult and of obtaining the M03 result is 



Piointi^k) ^qkB{M,N,pk) 



(;^)«'"»"-''- 

(19) 

At each redshift z^, we find the binomial prob- 
ability Pk that maximizes Pjoint(zfc). This set of 
optimal probabilities {pk} can be thought of as 
defining the parent population V of groups that 
maximizes the probability that we would obtain 
our null result and that the M03 groups would 
have their observed properties. Each percentage 
Pk in turn can be thought of as an estimate of the 
probability that a group would be visible in our 
survey if it were randomly selected from V and 
placed at redshift Zk- The value qk = ^ Pk), 
then, is the probability that such a group would 
not be detected in our survey. 

If 21 groups were selected at random from V 
and placed at the 21 redshifts Zk, the probabil- 
ity that none would be seen in our survey is the 
product of the 21 quantities qk- A maximal esti- 
mate of the conditional probability that none of 
the CN0C2 groups would be detected at the 5a 
level, given that the CN0C2 groups were drawn 
from the same population as the M03 groups, may 
therefore be represented as follows: 



21 



P{0 I same pop.) < J| (1 - Pfc) , 



(20) 



k=l 



where denotes our null result of not detecting 
any CN0C2 groups. The result of our simulation 
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is that the product in equation ([20)1 is P < 14.7% 
for M = 109 and P < 2.8% for M = 61; this indi- 
cates that the probabihty is small, but not vanish- 
ing, that the CN0C2 groups in our survey have 
the same luminosity function as the M03 groups. 

REFERENCES 

Akahori, T. & Masai, K. 2005, PASJ, 57, 419 

Akritas, M. G. & Bershady, M. A. 1996, ApJ, 470, 
706 

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 
2005, A&A, 441, 893 

Carlberg, R. G., Ycc, H. K. C., Morris, S. L., Lin, 
H., Hall, P. B., Patton, D. R., Sawicki, M., & 
Shepherd, C. W. 2001a, ApJ, 563, 736 

— . 2001b, ApJ, 552, 427 

Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 
28, 215 

Ettori, S., Tozzi, P., & Rosati, P. 2003, A&A, 398, 
879 

Fang, T., Gerke, B. F., Davis, D. S., Newman, 
J. A., Davis, M., Nandra, K., Laird, E. S., Koo, 
D. C., Coil, A. L., Cooper, M. C, Croton, D. J., 
& Yan, R. 2006, astro-ph/0608382 

Fukugita, M. & Peebles, P. J. E. 2004, ApJ, 616, 
643 

Giacconi, R., Murray, S., Gursky, H., Kellogg, E., 
Schreier, E., Matilsky, T., Koch, D., & Tanan- 
baum, H. 1974, ApJS, 27, 37 

Hoekstra, H., Franx, M., Kuijken, K., Carlberg, 
R. G., & Yee, H. K. C. 2003, MNRAS, 340, 609 

Hoekstra, H., Franx, M., Kuijken, K., Carlberg, 
R. G., Yee, H. K. C, Lin, H., Morris, S. L., Hall, 
P. B., Patton, D. R., Sawicki, M., & Wirth, 
G. D. 2001, ApJ, 548, L5 

Mahdavi, A. & Geller, M. J. 2004, ApJ, 607, 202 

McCammon, D., Almy, R., Apodaca, E., 
Bergmann Tiest, W., Cui, W., Deiker, S., 
Galeazzi, M., Juda, M., Lesser, A., Mihara, T., 
Morgenthaler, J. P., Sanders, W. T., Zhang, J., 
Figueroa-Feliciano, E., Kelley, R. L., Moseley, 
S. H., Mushotzky, R. F., Porter, F. S., Stable, 



C. K., & Szymkowiak, A. E. 2002, ApJ, 576, 
188 

McCammon, D. & Sanders, W. T. 1990, ARA&A, 
28, 657 

Mulchaey, J. S. 2000, ARA&A, 38, 289 

Mulchaey, J. S., Davis, D. S., Mushotzky, R. F., 
& Burstein, D. 2003, ApJS, 145, 39 

Osmond, J. P. F. & Ponman, T. J. 2004, MNRAS, 
350, 1511 

Ota, N. & Mitsuda, K. 2004, A&A, 428, 757 

Pearson, K. 1901, Philosophical Magazine, 2, 609 

Phonis, M. & Tovmassian, H. M. 2004, A&A, 416, 
441 

Press, W. H., Flanncry, B. P., Teukolsky, S. A., 
& Vetterhng, W. T. 1992, Numerical Recipes: 
The Art of Scientific Computing, 2nd edn. 
(Cambridge (UK) and New York: Cambridge 
University Press) 

Press, W. H. & Schechter, P. 1974, ApJ, 187, 425 

Rosati, P., Borgani, S., & Norman, C. 2002, 
ARA&A, 40, 539 

Rowan-Robinson, M. & Fabian, A. C. 1975, MN- 
RAS, 170, 199 

Scharf, C. 2002, ApJ, 572, 157 

Schwartz, D. A. 1978, ApJ, 220, 8 

Willis, J. P., Pacaud, F., Valtchanov, L, Pierre, 
M., Ponman, T., Read, A., Andreon, S., Altieri, 
B., Quintana, H., Dos Santos, S., Birkinshaw, 
M., Bremer, M., Due, P.-A., Galaz, G., Cosset, 
E., Jones, L., & Surdej, J. 2005, MNRAS, 363, 
675 

Xue, Y.-J. & Wu, X.-P. 2000, ApJ, 538, 65 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.0. 



16 



